Quintessence and phantom cosmology with non-minimal derivative coupling 
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We investigate cosmological scenarios with a non-minimal derivative coupling between the scalar 
field and the curvature, examining both the quintessence and the phantom cases in zero and con- 
stant potentials. In general, we find that the universe transits from one de Sitter solution to another, 
determined by the coupling parameter. Furthermore, according to the parameter choices and with- 
out the need for matter, we can obtain a Big Bang, an expanding universe with no beginning, a 
cosmological turnaround, an eternally contracting universe, a Big Crunch, a Big Rip avoidance and 
a cosmological bounce. This variety of behaviors reveals the capabilities of the present scenario. 



I. INTRODUCTION 

The cosmological research of the last three decades has 
elevated the role of scalar fields in the description of vari- 
ous sides of nature. Introduced as the driving mechanism 
for almost all inflation realizations the dynamics of 
scalar fields gained new interest after observations pro- 
vided indications for an accelerated universe expansion 
0. In particular, the new concept of "dark energy" was 
easier to be described by an extra scalar field dubbed 
quintessence 01 , than with the traditional cosmological 
constant [3, [ah an d the corresponding cosmological be- 
havior proves to be much richer. 

However, although the general belief is that the data 
are far from being conclusive, some data analyses sug- 
gested that the cosmological constant boundary, that is 
the phantom divide, has been crossed in the near cosmo- 
logical past [6|. The simplest way to explain this unex- 
pected behavior is the use of a phantom scalar field in- 
stead of a canonical one, that is a scalar with a negative 
sign of the kinetic term in the Lagrangian [7]. Although 
the discussion about the construction of quantum field 
theory of phantoms is still open in the literature (see for 
instance [8j for the causality and stability problems of 
phantom fields, but also [9[ for attempts in constructing 
a phantom theory consistent with the basic requirements 
of quantum field theory, with the phantom fields arising 
as an effective description), the richness and capabilities 
of phantom cosmology has gained significant interest in 
the literature, extended apart from dark energy to infla- 
tion area too 

Apart from the aforementioned basic use of scalars 
(canonical or phantom ones), cosmological models where 
the fields are non-minimally coupled to gravity [n], Il2j 
have been shown to present significant cosmological fea- 
tures, both for inflation and dark energy areas, and have 
been widely studied [l3j]. Additionally, one can further 



extend these "scalar-tensor" theories, allowing for non- 
minimal couplings between the derivatives of the scalar 
fields and the curvature [3], and these scenarios reveal 
interesting cosmological behaviors [lBj]. 1 

In our recent work we examined the cosmological 
scenario of a quintessence field with non-minimal deriva- 
tive coupling, and we extracted exact solutions in the case 
of zero potential. In the present work we are interested 
in extending this analysis in the case of non-zero poten- 
tials, and furthermore, for completeness, to perform it 
for both a quintessence and a phantom field. The plan of 
the manuscript is as follows: In section [TT] we construct 
the scenario and we extract the cosmological equations. 
In section Mil we examine specific potential choices and 
we investigate the corresponding cosmological solutions 
for various parameter choices. Finally, in section IIVI we 
discuss the physical implications of the different universe 
evolutions, and we summarize the obtained results. 



II. COSMOLOGY WITH NON-MINIMAL 
DERIVATIVE COUPLING 

In this section we present the cosmological paradigm 
with non-minimal derivative coupling between a scalar 
field and the curvature. In order to describe the 
quintessence and the phantom field in a unified way we 
adopt the e-notation, that is the parameter e takes the 
value +1 for the canonical field and —1 for the phantom 
one. 



A. Action and field equations 

Let us construct a gravitational theory of a scalar field 
with non-minimal derivative couplings to the curva- 



1 It is also worth mentioning a series of papers devoted to a non- 
*Electronic address: msaridak@phys.uoa.gr minimal modification of the Einstein- Yang-Mills-Higgs theory 
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ture. In general one could have various forms of such 
couplings. For instance in the case of four derivatives one 



could have the terms KxRcf)^'^ , ^RuuC, 



K4,R[ iV <t>4>'V' v , K^R-^(j)(j>'^ and kqORc/) 2 , where the coeffi- 
cients ki,...,kq are coupling parameters with dimen- 
sions of length-squared. However, as it was discussed 
in [13, EH, El] , using total divergencies and without loss 
of generality one can keep only the first two terms. Thus, 
the action for the cosmological scenarios at hand writes: 



-9 



R 
8^ 



[eg^ + KG^] cj>^<j>>» -2V(<j>) 



(1) 

where V{4>) is a scalar field potential, g^ is a metric, g — 
det(g^), R is the scalar curvature, G M „ is the Einstein 
tensor, and k is the single derivative coupling parameter 
with dimensions of length-squared. 

Varying the action ([1]) with respect to the metric g^ 
leads to the gravitational field equations 



Gfiu = 8n[eT^ + kQ„ v ] - torg^V (4>), (2) 



with 



Tfnu = V^V„^- ig^(V0) 2 , 

= -^V„<l>V v 4>R + 2V a <l>V i ^R^ 

-v M v^n</>-i(v^) 2 G^ 

+9»*[ - |V a V^V a V^ + i(D0) 2 

Similarly, variation of the action (JlJ with respect to <f> 
provides the scalar field equation of motion: 

[egT + kG^V^V^ = V+, (3) 
where = dV(<p)/d(p. 

B. Cosmological equations 

Throughout this work we consider a spatially-flat back- 
ground geometry with a metric 



ds 2 



-dt 2 + e 2a Mdx 2 , 



(4) 



where a(t) = e°''' is the scale factor, and dx 2 is the 
Euclidian metric. Thus, the Hubble parameter is simply 
H(t) = a{t)/a{t) = d(i). 

As usual we assume a homogenous scalar field, namely 
4> = 4>{t). In this case the field equations ([2]) and ([3]) are 
reduced to the following system: 



3d 2 = 4tt</> 2 (e 



9na 



-2a. - 3d 2 = 4tt0 2 



+ 8ttV(^>) 
k (2a 



3d 2 



(5) 



■8ttV(4>), (6) 
V*, (7) 



e((f> + 3d<j>) - 3K,(a 2 (f> + 2dd0 + 3d 3 </> 
where a dot denotes a derivative with respect to time. 
Note that equations ((6]) and © are of second order, while 
([5]) is a first-order differential constraint for a(t) and <fr(i). 

The constraint ([5]) can be rewritten as: 



3d 2 - 87rV(0) 
4tt(£ - 9na 2 ) ' 



(8) 



or equivalently as 



4ire<p 2 + 8ttV((/)) 

3(1 + 127TK(/) 2 ) 



(9) 



Therefore, as long as the parameters e and re and the 
potential T^(<^) are given, the above relations provide re- 
strictions for the possible values of a and <f>, since they 
have to give rise to non- negative 4> 2 and a 2 , respectively. 

Let us now separate the equation for (j) and a. For this 
aim we resolve equations ^ and (J7J with respect to a 
and </) and, using the relations (jHJ and (|9|), we eliminate 
and d from respective equations. We easily result to: 



-2V3^<j)[e + s8iTK<j) 2 - 8TrnV{(j))]^[s<i) 2 + 2^0)](127rre ( /> 2 + 1) + (12ttk0 2 + l)(47rre</> 2 + 

£(1 + 127TK0 2 + 967T 2 K 2 (/) 4 ) + 87rKy(0)(127TK(/) 2 - 1) 

.. _ -(e - 3Kd 2 )(e - 9«d 2 )[3d 2 - 8ttV(0)] - ±ifxK,a^{e - 9«;d 2 )[3d 2 - 8irV(<p)} V4, 
1 — 9ered 2 -I- 54« 2 d 4 — 8nKV(<j))(e + 9na 2 ) 

I 



We mention however that although the ^-equation does not contains a-terms, the a-equation in general contains 
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-terms arising from the potential V(4>). 



III. COSMOLOGICAL SCENARIOS AND 
SOLUTIONS 

In this section we examine specific cosmological sce- 
narios, that is we consider specific potential choices. 



A. Zero potential: V(<j>) = 

The case of a canonical field under zero potential has 
been investigated in [3]. Thus, in this subsection we 
restrict ourselves in the case of a phantom field, that is 
e = — 1. The field equations ©-([7]) now read 



3d 2 



-4tt</> 2 (l + 9/td 2 ) 



-2d - 3d 2 = -4tt</> 2 



1 - k 2d + 3d 2 



3ad> + 3k 



a 2 4> + 2aa<j) + 3d 3 </> = 

The constraint (fT2|) can be rewritten as: 

3d 2 



0. 



(12) 

4d0<T 1 ) (13) 
(14) 



4vr(l + 9Kd 2 ) : 



or equivalently as 



47T0 2 



3(1 + Uirnq 



(15) 



(16) 



From these relations we deduce that d and <f> should obey 
the following conditions: 



1 + 9«:d 2 < 

1 + 127TK</> 2 < 0. 



(17) 
(18) 



Note that these conditions are only fulfilled for n < 0. 
Assuming k = — k , wc find 



d 2 > — 



> 



(19) 
(20) 



The separate 4> and a-equations f ([T0|) and (|TT|) respec- 
tively) read: 



2V37r0 2 (87rK<? 



1)\ 12ttk6 2 - 1 



1 - 12ttk0 2 + 96tt 2 k 2 4 



3d 2 (3fc 2 d 2 - l)(9fc 2 d 2 - 1) 
1 - 9fc 2 d 2 + 54fc 4 d 4 



(21) 



(22) 



From these equations we deduce that the case at hand 
exhibits three qualitatively different cases: 

Al. (9k 2 )^ 1 < a 2 < (3/c 2 ) -1 . In this case d is posi- 
tive, and d increases with time. The solution a is varying 



between two de Sitter asymptotics: at_ i ._ 00 = t/3k and 
a t ^oo = t/V3k. 

A2. d 2 = (3fc 2 ) -1 , d = 0. One has exactly the de 
Sitter solution 

a(t) = t/V3k, <j)(t) = t/Vs^k. 

A3, d 2 > (Sfc 2 )" 1 . In this case d is negative, and 
d decreases with time. Nevertheless, we remind that d 
remains larger than 1/^/Sk always. The solution a(t) is 
varying between two asymptotics. The t — > oo asymp- 
totic is the de Sitter one: a t ^oo — t/v3k. The second 
asymptotic can be obtained as follows: Assuming d — > oo 
at t —> ti gives the following asymptotical form of Eq. 
(123): 



—a 

2 



12tt/c 2 



with the asymptotic solution 



= ai + - ln(t - U) 



2V3irk 



(t-U). 



(23) 
(24) 

(25) 
(26) 



Thus, by construction, the moment t = ti corresponds to 
a cosmological singularity. 

In order to present this cosmological behavior more 
transparently, we perform a numerical elaboration of the 
model at hand, namely of equation (|2"2"j) . and the results 
are presented in Fig. [TJ We mention that the possible 




(a) 



(b) 



FIG. 1: The evolution of a(t) in the phantom case with zero 
potential. The dash-dotted and dotted lines denote the de 
Sitter asymptotics t/3k and t/y3k } respectively. The solid 
lines corresponds to a(t) with (a) (9A: 2 ) -1 < d 2 < (3/c 2 ) -1 
(case Al); (b) a 2 > (3fc 2 ) -1 (case A3). We have considered 
k 2 - -L 

— 27- 

singular behavior of a(t) means that the scale factor itself 
(a(t) = e a ^) starts from zero at some initial time. 
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B. Cosmological constant: V((j>) = A = const 



In this particular scenario the separate 
equations (ITU1) and (ITTj) read: 



and 



8ttkA] 



127TK(/) 2 + 96tT 2 K 2 4 ) + 87TKA(127rK0 2 - 1) 

(27) 



3(e A Hl - d 2 )(e - 3Kd 2 )(e - 9na 2 ) 
1 - 9eKd 2 + 54k 2 « 4 - 3ne A Hl(e + 9na 2 ) ' 



(28) 



where for simplicity we have defined -ff a = -y/ 8vr | A. | / 3 and 
e A = sign A. It proves convenient to consider separately 
the various cases arising from specific choices of the pa- 
rameters e, A, and k, and in the following we present the 
eight qualitatively different cases of the scenario at hand. 

Bl. e = 1, A> 0, k > 0. 

It is easy to see that the a-equation |28|) has three 
trivial particular solutions: (i) a(t) — H A t, (ii) a(t) = 
t/y/SK, and (iii) a(t) = t/y/9~K. Sequentially, substituting 
them into the whole system of field equations (HJ)-®, one 
may straightforwardly find the following exact solutions: 



a(t) 
a(t) 



H A t, 
t 



0(t) = 



4>q = const, 



1 



8lTK 



t, H A > 



'3k 



which describe de Sitter universes. However, the third 
solution a(t) = t/\9~K cannot in general satisfy equations 
©-© (apart from the fine-tuned case l/y9~K = H A ). 

More generally, the constraint ((8]) gives the following 
restrictions for d 2 : 



x\ < a < X2, 



(29) 



where X\ = min(l/9/t, H^) and X2 = max(l/9«;, ifjQ 



The second derivative d, given by relation (|28p . is nega- 
tive if H\ < 1/9k, and positive if H\ > 1/9k, and thus 
d is decreasing or increasing with time. As it is deduced, 
the corresponding solutions for a(t) are varying between 
the two de Sitter asymptotics depending on values of pa- 
rameters H A and n. This behavior can be observed in 
Fig. [2j arisen from numerical elaboration. 

B2. £ = 1, A> 0, k < 0. 

In this case equations ©-© have the de Sitter solu- 
tion a(t) = H A t, <p(t) = 4>o = const. Additionally, the 
constraint (HJ) yields 

a 2 > H%. 

Under this condition d, given by relation (|28[). is nega- 
tive and d is decreasing with time. The corresponding 
solution a(t) is singular at some initial moment of time, 
i.e., lim t _,. ti a(t) = — oo, while for large times a(t) tends 
to de Sitter asymptotic H A t. This behavior can be seen 
in Fig. H 
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(c) 



FIG. 2: The evolution for a(t) for e = 1 (quintessence), 
A > 0, k > (case Bl). The dashed, dash-dotted, and dotted 
lines denote de Sitter asymptotics a(t) = HAt, t/V9/t, and 
t/v3ft, respectively. A solid line corresponds to a(t) with (a) 
Hi < a 2 < 1/9k; (b) 1/9/t < a 2 < Hi < 1/3/s; (c) 1/9k < 
d 2 < 1/3/t < H%. In graph (a) we have considered Hi = 1, 
in graph (b) Hi = 5, and in graph (c) Hi = 10; everywhere 



B3. £ = -1, A> 0, k> 0. 

For these parameter choices, equations ©-([7]) have the 
de Sitter solution a(t) = H A t, <p(t) = 4>o — const. Gen- 
erally, the constraint © gives 

a 2 < H 2 . 

Under this restriction d, from (|28[) . is positive and d is 
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FIG. 3: The evolution for a(t) for s = 1 (quintessence), 
A > 0, K < (case Bl!,). The dashed line corresponds to 
the de Sitter solution a(t) = H^t with Ha < 1/ \/3|re| . The 
solid curve corresponds to a solution with a 2 > H A . We have 
considered H\ = 1 and k ■ 



j_ 

' 27 ' 



increasing with time. The corresponding solution for a(t) 
varies between two de Sitter asymptotics: a(t — > — oo) ~ 
—H\t and a(t — ¥ oo) w H\t. This behavior is more 
transparently shown in Fig. |4j 




Ha and k. The corresponding possible types of solutions 
are demonstrated in Fig. [5] 
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(e) 



FIG. 4: TTie evolution for a(t) for e = — 1 (phantom), A > 
0, k > ('case TTie dashed line corresponds to the de 

Sitter solution a(t) = Htd with Ha < 1/\/~3k. The solid line 
corresponds to a solution with a 2 < H\. We have considered 
Ha = 1 and k = -k~ ■ 



B4. e = -1, A> 0, k < 0. 

In this case equations ©-(O possess two different de 
Sitter solutions: 



const, 




'3|k|' ' " V 8?r N 
The constraint ([8]) now yields 



3 ^t, m < 



a < x\ 



a > x 2 , 



where x\ — min(l/9/-c, H\) and x 2 — max(l/9ft, H%). 
Thus, the behavior of a(t) satisfying the above condi- 
tions depends on the specific values of the parameters 



FIG. 5: The evolution of a(t) for e = — 1 (phantom), A > 0, 
ft < (case B4)- The dashed, dash-dotted, and dotted lines 
denote the de Sitter asymptotics a(t) — HAt, tj \J9\k\ , and 
tj \J?>\k\, respectively. The solid lines corresponds to a(t) with 
(a) Hi < 1/9H < cV 2 < 1/3|k|; (b) Hi < 1/9\k\ < 1/3\k\ < 



(c) l/3k| 
< Hi 



<H 2 



< a 2 : 



(d) a 2 < H 2 A < 1/9\k\; (e) a 2 < 



1/9|«| 

H\ — 1, in graphs (c) and (e) H\ = 15; everywhere k = — 



In graphs (a), (b) and (d) we have considered 

— i . t " ' ' ; 1 i- i i i t i ■ / 1 i ■ ri ;-: — — 

27 ' 



B5. e = 1, A< 0, k > 0. 

In this case the constraint © leads to 



d 2 < 



1 

9k' 



Under this condition a, from relation (|28|). is negative 
and a is decreasing with time. The corresponding so- 
lution for a(t) varies between the two de Sitter asymp- 
totics: ctt->-oo ~ i/v9K and a^oo ~ — i/v9K- These 
features are presented in Fig. [51 arisen from numerical 
elaboration. 
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FIG. 6: The evolution of ct(t) fore = 1 (quintessence) , A < 0, 
K > (case B5). The dash-dotted line denotes the de Sitter 
asymptotic a(i) = t/y§K. The solid curve corresponds to a(t) 
with a 2 < 1/9k. We have considered H\ — 1 and k = 



B6. e = 1, A< 0, k < 0. 

For these parameter sub-class, the constraint ([8]) does 
not lead to any restriction on the values of a 2 . Addi- 
tionally, the second derivative a, from ([25]). is negative 
and a is decreasing with time. The corresponding solu- 
tion for a(t) varies between the two singular solutions: 
a t ^±f = — oo. This behavior can be seen more trans- 
parently in Fig. [7J 



The constraint ([8]) now gives 



Under this condition a from (|28p is positive if l/9|rt| < 
a 2 < 1/3 1 k\, and negative if a 2 > l/3|re|; respectively, a 
is increasing or decreasing with time. The two possible 
types of solutions for a(t) are presented in Fig. |SJ arisen 
from numerical elaboration. 
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(a) (b) 

FIG. 8: The evolution of a(t) for e = —1 (phantom), A < 0, 
k < (case B8). The dash-dotted, and dotted lines denote de 
Sitter asymptotics a(t) — t/W9\n\, and t/^/3|/t|, respectively. 
The solid curve corresponds to a(t) with (a) 1/9|k| < a 2 < 
l/3|/s|; (b) a 2 > l/3|/t|. We have considered H\ — 1 and 
K = - 27 ■ 
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We close this subsection by mentioning that for A — > 0, 
the scenario at hand coincides with the previously studied 
quintessence case with V = [HI], or with the phantom 
with V = examined in the previous subsection. 

Finally, note that in principle one could extend the 
aforementioned analysis to more complicated potentials, 
such as the exponential and the power-law one. Unfor- 
tunately, the complexity of the non-minimal derivative 
coupling does not allow for the extraction of any ana- 
lytic solutions in these cases. The examination of the 
cosmological behavior of such scenarios must be based 
on numerical investigation, and this is left for a future 
investigation [l9j . 



FIG. 7: The evolution of a(t) for e = 1, A < 0, k < (case 
B6). We have considered H A — 1 and k = — 



IV. DISCUSSION AND CONCLUSIONS 



B7. £ = -1, A< 0, k > 0. 

In this particular case, the constraint © cannot be 
fulfilled. Therefore, there are no solutions corresponding 
to this scenario sub-class. 

B8. £ = -1, A< 0, K < 0. 

In this case equations ©-([7]) posses the following de 
Sitter solution: 



*(*) = 



t 




3MH 2 



8tt« 



(30) 



In this work we investigated cosmological scenarios 
where there is a non-minimal derivative coupling be- 
tween the scalar field and the curvature. In order to be 
complete, we considered both quintessence and phantom 
fields, although the later case could be ambiguous at the 
quantum level. Finally, in order to examine the pure ef- 
fects of these scenarios, we have not included the matter 
content of the universe, although this can be straightfor- 
wardly taken into account. 

A first observation is that the non-minimal derivative 
coupling leads to qualitatively different behavior, com- 
paring to the uncoupled case, even for the simple cases 
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of zero or constant potentials. In particular, as we ob- 
serve the universe evolves between two asymptotic de 
Sitter solutions, characterized by the strength of the cou- 
pling. In the limit where the coupling tends to zero and 
in the cases where the solutions exist, these two asym- 
potics coincide and the system acquires only one de Sitter 
solution, namely the one that is exhibited from the cor- 
responding conventional (that is uncoupled) model of a 
universe without matter. We mention that the transi- 
tion between the two de Sitter solutions is a pure effect 
of the non-minimal derivative coupling, and it does not 
require the presence and the role of matter. Finally, note 
that as time roll backwards, the scale factor of the uni- 
verse (a(t) — e a (*)) can be either eternally decreasing, or 
become zero at some initial time. Thus, our scenario 
exhibits either the Big Bang, or it corresponds to an 
eternally expanding universe with no beginning, with the 
later case arising easily, without the need of a specially 
designed potential as in conventional cosmology 20] . 

An additional feature of the scenario at hand is the 
radically different evolution of a quintessence universe in 
some solution sub-classes. In particular, for negative cos- 
mological constant and positive coupling (case B5) the 
scale factor of the universe is growing, it reaches a maxi- 
mum, and then it decreases. This is the realization of the 
cosmological turnaround, in which the universe transits 
from expansion to contraction 0, [2l| . The fact that this 
is obtained solely from the dynamics of the non-minimal 
derivative coupling, without the need for matter or for 
exotic gravitational terms, makes the scenario at hand 
very interesting. Lastly, note that the contracting phase 
is eternal, that is the universe does not result to a Big 
Crunch [|. 

In similar lines, in the quintessence case with negative 
cosmological constant and negative coupling (case B6), 
the scale factor starts from zero at some initial time and 
returns to zero at some final time. This is the realization 
of a universe starting from a Big Bang and ending with a 
Big Crunch [5| , and the fact that this is obtained without 
the need of matter is a novel effect of the non-minimal 
derivative coupling. 

The phantom evolution exhibits also the aforemen- 
tioned behavior. Apart from an eternally expanding uni- 
verse, including the transition between two de Sitter so- 
lutions, as we observe in Fig. 5(e) (case B4), that is 
for positive cosmological constant and negative coupling, 



the universe can experience the cosmological turnaround. 
This is radically different comparing to the uncoupled 
phantom scenarios, which not only cannot experience the 
turnaround, but on the contrary, in the absence of mat- 
ter they result to a Big Rip [22]. It seems that the non- 
minimal derivative coupling smoothers or (for large cou- 
pling) completely alters the evolution, leaving phantom 
cosmology free of a Big Rip. This new and significant 
behavior reveals the richness of the scenario at hand. 

However, the phantom case exhibits an additional sur- 
prising feature, that is not present in the quintessence 
scenario. As we observe in Fig. 4 (case B3, that is posi- 
tive cosmological constant with positive coupling) , as well 
as in Fig. 5(d) (case B4, that is positive cosmological 
constant with negative coupling), the universe can tran- 
sit from the contracting to the expanding phase, without 
meeting any singularity. This is just the cosmological 
bounce [23j , and its realization from a sole phantom field 
make the scenario at hand very interesting. 

In summary, the paradigm of non-minimal derivative 
coupling either in the quintessence or in the phantom 
case, may have important cosmological implications, even 
in its simplified realization where matter is absent. Apart 
from the transition between different de Sitter solutions, 
according to the parameter choices we can obtain a Big 
Bang, an expanding universe with no beginning, a cos- 
mological turnaround, an eternally contracting universe, 
a Big Crunch, a Big Rip avoidance and a cosmological 
bounce, and this variety of behaviors reveals the capabil- 
ities of the scenario. Furthermore, one could generalize 
this paradigm in the case where both the quintessence 
and the phantom fields are present, that is to generalize 
the so called "quintom" paradigm [24j in the case of non- 
minimal derivative coupling. In these scenarios one could 
obtain the combination of the above behaviors, such as 
to obtain a cyclic cosmology [25[. Definitely, this subject 
deserves further investigation. 
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